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Radiation force balance (RFB) measure- 
ments of time-averaged, spatially-integrat- 
ed ultrasoimd power transmitted into a 
reflectionless water load are based on 
measurements of the power received by 
the RFB target. When conical targets are 
used to intercept the output of coUimated, 
circularly symmetric ultrasound sources 
operating at frequencies above a few 
megahertz, the correction for in-situ atten- 
uation is significant, and differs signifi- 
cantly from predictions for idealized cir- 
cumstances. Empirical attenuation correc- 
tion factors for a 45° (half-angle) absorp- 
tive conical RFB target have been deter- 
mined for 24 frequencies covering the 5 
MHz to 30 MHz range. They agree well 
with previously impubhshed attenuation 
calibration factors determined in 1994 for 
a similar target. 
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1. Introduction 

Radiation pressure [1-13] has been employed in a 
wide variety [14-27] of ultrasound power meter 
designs. In these instruments, an appropriately con- 
structed target properly aligned in a steady-state under- 
water ultrasound field is subjected to a radiation force 
F given by 



F=erWlc 



(1) 



where Cj. is a correction factor determined by various 
properties of the target, W is the time-averaged spatial- 
ly integrated power intercepted by the target, and c is 
the speed of sound in the water Time-averaging of the 
ultrasound occurs because, under practical circum- 



stances, the inertia of the target causes it to effectively 
integrate pulses into the corresponding steady-state 
force. Spatial integration is a consequence of the 
extended geometry of all targets and is optimized in 
practical instruments by using targets larger in cross 
section than the incident beam. 

One such instrument was designed and custom built 
at the National Bureau of Standards (NBS) to be a pri- 
mary standard for ultrasound power, and served this 
purpose for the 24 years (1977-2001) that NBS/NIST 
ultrasound power measurement services were available 
to the public. The features which distinguish its design 
from all other known RFB designs, and the evolution of 
the means used to implement these features, have been 
described extensively elsewhere [23,26-27]. This arti- 
cle presents a detailed description of the method used 
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since 1994 to correct for the attenuation of ultrasound 
within the NIST RFB in the configuration used to 
measure the output of ultrasound power transfer stan- 
dards [29-30]. Empirical attenuation correction factors 
determined in 2005 are presented for 24 frequencies 
covering the 5 MHz to 30 MHz range, and are com- 
pared with previously unpublished attenuation calibra- 
tion factors determined in 1994. 

2.1 RFB Targets 

Selection of the optimal shape and material for an 
RFB target depends on the usual variety of application- 
specific details. Combinations of constraints on accura- 
cy, cost, propagation medium, ultrasound amplitude, 
frequency, spatial (scanned v. stationary), and temporal 
(continuous v. pulsed) characteristics have led to the 
use of flat, waffled, concave, and convex shapes, and to 
the use of both absorptive and reflective materials [14- 
27]. 

A useful feature of the conical target shape was iden- 
tified by Borgnis [5], who found that for a conical tar- 
get of half-angle (j) centered in an unfocused ultrasound 
beam with propagation axis parallel to that of the cone, 



ej-= I -P^ cos If 



(2) 



where j3 is the amplitude reflection coefficient of the 
target. From Eq. (2) it follows that ej- can be set to unity 
for all values of p by the simple expedient of setting (j) 
to 45°. 

Results of experimental research begun in 1974 at 
NBS have confirmed that gj- can be made to differ from 
unity by less than 0.1 % by setting (j) to 45° for targets 
made of metal or absorptive rubber. Because targets of 
this type were the only ones used in the work described 
in this article, no further mention of gj- will be made. 

2.2 In-situ Attenuation with tlie 45° Conical 
Target 

Most applications require that RFB measurements be 
used to determine the total power Wq radiated from the 
output port of a transducer, rather than the power W 
intercepted by the RFB target. 

For the RFB configuration depicted in Fig. 1, let P be 
the average ultrasound pressure in the plane which 
intersects the tip of the target and is perpendicular to the 
axis of the ultrasound beam. With a collimated, circu- 
larly symmetric ultrasound source of radius R located a 
distance d from this plane. 



FORCE SENSOR 




ULTRASOUND 
SOURCE 



Fig. 1. Simplified configuration of NIST RFB. Shading indicates sil- 
icone rubber conical shell of target. 



p^Ke- 



cnR^ 



(3) 



where a is the empirically determined attenuation coef- 
ficient for the distilled water propagation medium, and 
c is the speed of sound in distilled water. An implicit 
condition for the validity of Eq. (3) is that a be inde- 
pendent of position in the ultrasound beam between the 
transducer and the target. This has long since been ver- 
ified by experimental data, and is confirmed independ- 
ently by recent data presented later in this article. 

For the ultrasound sources used in this study, P is 
assumed to be uniform within the ultrasound beam. 
Because (/>=45°, the perpendicular distance from the 
apex plane to a point on an annulus on the target sur- 
face is identical to the radius of that annulus. Taking all 
annuli into account, the radiation force F induced on 
the target is given by 



F= {" Pe-^'lnrdr 

JO 



which reduces to 



F=^\l-(aR+\)e''n 



(4) 



(5) 



436 



Volume 111, Number 6, November-December 2006 

Journal of Research of the National Institute of Standards and Technology 



Combining Eqs. (1), (3), and (5), we define attenuation 
correction factor C as 



c^w^/w^ 



a'R' 



2e-"''[l-e-'"'(aR + l)] 



(6) 



3. Experimental Procedure 

In-situ attenuation was determined empirically by 
measuring the output of a NIST Standard Ultrasound 
Source (SUS) [29-30] for 1 1 values of d covering the 
range 1 mm to 11 mm in 1 mm increments indicated 
with approximately 0.1 mm resolution. The target for 
all tests was a 45° half- angle conical shell of silicone 
rubber, 10 mm thick (measured axially) and 52 mm in 
diameter, surrounding an inner cone of rigid 
polyurethane foam. An ultrasound ray arriving parallel 
to the axis of the cone is reflected only slightly at the 
water-rubber interface because the rubber offers a good 
acoustical impedance match to that of water In passing 
through the rubber, the ultrasound loses most of its 
energy to absorption. Near total reflection occurs at the 
interface between the rubber and the polyurethane 
foam, and the reflected ultrasound is further attenuated 
as it travels outward radially through the rubber layer 
With most of the incident ultrasound absorbed, the 
remainder is reflected radially outward toward the 
cylindrical wall of the test tank. This wall is lined with 
irregularly oriented wedges made of AptFlex F36' 
acoustic absorbing rubber These wedges are approxi- 
mately 20 mm by 25 mm at the base, approximately 35 
mm long, and extend radially inward from the cylindri- 
cal wall. 

The SUS transducer was installed in a customized 
test tank fitting which allows direct coupling to the 
water in the test tank with no intervening membrane. 
The transducer was positioned so that its output plane 
was several millimeters above the bottom surface of the 
tank. This arrangement ensures that any radiated edge 
waves will impinge on the absorptive wedges instead of 
being scattered from the bottom surface of the tank. 

For all tests, the radio frequency drive voltage was 
adjusted to obtain 5.000 volts at the SUS dc monitor 
output [30]. This corresponds to output power levels 
equivalent to roughly 40 % of the rated SUS maximum 



Certain commercial equipment, instruments, or materials are iden- 
tified in this paper to foster understanding. Such identification does 
not imply recommendation or endorsement by the National Institute 
of Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the pur- 
pose. 



power output [29-30]. Variations in drive voltage 
greater than 0.02 % were prevented by custom-built 
automatic-level-control circuitry [26]. Because the 
radiation conductance of the SUS transducer varies 
with frequency, the transducer power output varied 
from 348 mW to 455 mW for the 4.8 MHz to 29.2 MHz 
test frequency range. The test frequencies correspond- 
ed to odd harmonics of the approximately 0.5 MHz 
fundamental resonance frequency of the SUS transduc- 
er, and covered the test frequency range in increments 
of approximately 1 MHz. 

For each of 1 1 values of liftoff distance d, the water 
temperature Twas measured and recorded and the SUS 
dc monitor voltage was verified. Then 5 measurements 
of the intercepted power W were made. Each measure- 
ment was made by manually adjusting the magnitude 
and phase of a nulling signal to achieve an adequate 
null as indicated by an oscilloscope and galvanometer 
The nulling voltage was then recorded in a spreadsheet 
which used the test temperature and an RFB scale fac- 
tor to calculate the corresponding value of W. 

System stability was monitored for each set of 5 
measurements by calculating o/|i, the relative standard 
deviation of the mean. The lower limit of this parame- 
ter for NIST RFB measurements made under typical 
circumstances [27] is 0.8 %. By this measure, system 
performance exceeded expectation during the course of 
the work reported in this article, as the largest value of 
o/|i for any data set was 0.31 %. 

The time of day was logged at the beginning of each 
set of measurements. The average time required to 
make each of 1320 measurements was 29 seconds. 
Results for a typical test are shown in Table 1 . 



4. Data Analysis 

For each set of 1 1 power measurements for a partic- 
ular ultrasound frequency, the in-situ attenuation coef- 
ficient of water was calculated from the slope of a lin- 
ear regression fit oi\n{W„^JW^ to d, where W^ is the 
average power intercepted by the RFB target with d at 
its minimum value. Results for the data of Table 1 are 
plotted in Fig. 2. The linear least-squares regression 
trendline in the figure lies within the bounds of error 
bars that reflect the 2 % best-case combined relative 
expanded power measurement uncertainty [27]. 

The results for all test frequencies are shown in Table 
2. Attenuation coefficient a, in nepers per meter, is the 
slope of the trendline of the linear least-squares regres- 
sion fit oi\n{W„^JW^ to d, with d expressed in meters. 
The parameter R is the Pearson product-moment corre- 
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Table 1. Data for 17.495 MHz test conducted on 6 July 2005 



Time 


^("C) 


liftoff rf 


fVl 


W2 


fV3 


fV4 


W5 


W 

" avg 


EDT 




(mm) 


(mW) 


(mW) 


(mW) 


(mW) 


(mW) 


(mW) 


1341 


21.7 


1 


396.5 


397.0 


396.8 


396.7 


397.2 


396.8 


1344 


21.7 


2 


392.9 


393.2 


394.0 


393.7 


393.5 


393.5 


1349 


21.7 


3 


392.8 


393.5 


393.4 


391.6 


391.3 


392.5 


1353 


21.7 


4 


386.6 


386.6 


386.5 


386.5 


386.5 


386.5 


1355 


21.7 


5 


384.9 


385.0 


384.7 


384.7 


384.7 


384.8 


1356 


21.7 


6 


378.4 


379.5 


379.1 


379.2 


379.1 


379.1 


1358 


21.7 


7 


375.4 


375.3 


375.4 


375.4 


375.6 


375.4 


1359 


21.7 


8 


372.3 


372.6 


372.8 


372.8 


372.4 


372.6 


1401 


21.7 


9 


364.7 


364.7 


364.6 


364.2 


364.6 


364.6 


1402 


21.7 


10 


363.6 


363.5 


363.3 


363.3 


363.0 


363.3 


1404 


21.7 


11 


356.8 


357.2 


356.4 


357.2 


357.2 


357.0 
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Fig. 2. Variation of log-normalized power with liftoff distance,/=17.495 MHz. 



lation coefficient [28]. To make it easier to see how 
closely R approaches unity in all cases, the parameter 
\-R is tabulated instead. 

The variation of \-R with /in Table 2 indicates that 
the goodness of fit does not vary significantly over the 
range of/ The high values of i? for all test frequencies 
confirm that a was constant throughout the volume of 
water defined by the ultrasound beam between the 
transducer and the surface of the target. Because this 
volume changed significantly as liftoff distance d was 
increased by 10 mm, the correlation indicated by R 
would have been compromised if a had changed signif- 
icantly within this volume. 



4.1 Long-term Reproducibility 

The tests just described were conducted in June and 
July of 2005. To explore the effects of slight variations 
in system characteristics, we next examine the results 
of similar tests conducted in September and October of 
1994 using a target which was identical except that its 
inner cone was constructed of polystyrene foam, rather 
than the polyurethane foam used in the target tested in 
2005. The other systematic differences applicable in 
1994 were the use of a test tank equipped with wedges 
made of organic rubber, instead of AptFlex F36, and the 
use of four SUS transducers, rather than one. The 
results of these tests are presented in Table 3. 
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Table 2. NIST RFB in-situ attenuation coefficients, from tests con- 
ducted in 2005 

/(MHz) a(Np/m) \-R /(MHz) a(Np/m) \-R 



4.768 

5.828 

6.890 

7.952 

9.013 

10.072 

11.133 

12.193 

13.256 

14.315 

15.375 

16.436 



1.10 
1.37 
1.70 
2.63 
2.89 
4.22 
4.40 
5.18 
5.95 
6.85 
8.22 
9.27 



0.004 
0.011 
0.002 
0.004 
0.001 
0.010 
0.001 
0.003 
0.001 
0.002 
0.001 
0.001 



17.495 
18.556 
19.617 
20.679 
21.740 
22.799 
23.860 
24.920 
25.980 
27.041 
28.101 
29.161 



10.6 
12.1 
13.1 
14.7 
15.7 
17.4 
18.8 
20.7 
23.2 
24.2 
26.3 
28.9 



0.008 
0.001 
0.001 
0.001 
0.001 
0.002 
0.001 
0.002 
0.001 
0.001 
0.002 
0.001 



Table 3. NIST RFB in-situ attenuation coefficients, from tests con- 
ducted in 1994 



/(MHz) 


a (Np/m) 


\-R 


/(MHz) 


a (Np/m) 


\-R 


5.828 


1.32 


0.017 


20.679 


13.6 


0.003 


9.060 


2.40 


0.006 


21.140 


15.3 


0.001 


9.447 


2.76 


0.003 


22.248 


15.0 


0.002 


10.072 


3.60 


0.001 


24.980 


21.2 


0.001 


15.100 


7.41 


0.007 


25.980 


21.6 


0.001 


15.375 


8.29 


0.002 


27.205 


24.9 


0.001 


15.860 


8.77 


0.001 


30.220 


28.9 


0.002 



Table 4. Least-squares fit results for NIST RFB in-situ attenuation 

—1 —2 

coefficients. Slopes are expressed in Np m (MHz) , intercepts in 
Npm"' 



Date 

2005 
1994 



a to/ slope 

0.0331 
0.0322 



a to/ intercept 

0.338 
0.207 



1.000 
0.998 



ing wave motion that occurs when some of the ultra- 
sound energy lost to attenuation is translated into 
momentum, in effect inducing a radiation force on the 
water itself Because the NIST RFB uses audio-fre- 
quency pulse modulation of the ultrasound being meas- 
ured, the ultrasound energy lost to streaming takes the 
form of water waves from which the target cannot be 
isolated because isolating membranes with acceptable 
ultrasound transmission characteristics inevitably 
would offer negligible resistance to the transmission of 
audio frequency water waves. With the attenuation in 
water virtually nil for audio frequencies, the amplitude 
of these water waves is undiminished as they travel to 
the surface of the target to induce forces which add to 
those induced by the ultrasound itself, reducing the 
apparent attenuation of the ultrasound from the source 
[26]. 

4.3 Attenuation Correction Factor Results 



A cursory inspection suffices to establish that the 
variation of \-R with/in 1994 was about the same as it 
was in 2005. The calculated average values of \-R, 
0.0034 for the 1994 data and 0.0027 for the 2005 data, 
are insignificantly different. 

4.2 Frequency Dependence 

With a now available for a number of ultrasound fre- 
quencies determined by transducer resonances, we next 
consider interpolation schemes for calculating a for 
arbitrary frequencies within the 5 MHz to 30 MHz 
range. Methods involving fitting atof^ are attractive 
because the calculated slope is directly comparable 
with the published attenuation coefficient for distilled 
water [31]. Least-squares linear regression results for 
all data in Tables 2 and 3 are shown in Table 4. 

While only slightly different from each other, the 
two slopes are substantially smaller than 0.0472 Np m"' 
(MHz)~^, the attenuation coefficient for water at 22° C 
derived from values reported by Bacon and Jarvis [31] 
for distilled water and ideal (free-field plane wave) 
conditions. This discrepancy is attributed to the stream- 



We begin by examining the accuracy of the least- 
squares-fit interpolation method for determining values 
of a for arbitrary frequencies. Hereinafter, all values of 
C will be calculated for liftoff distance d=l mm, the 
value typical of NIST RFB operations, and fori?=7.938 
mm, the radius of the NIST SUS transducer 

For each of the 24 frequencies of the tests conducted 
in 2005, Fig. 3 shows values of C based on two values 
of a, one determined from the test results for that fre- 
quency, and the other determined by interpolation using 
the regression results of Table 4. For clarity, the bounds 
of the applicable 2 % combined relative expanded 
power measurement uncertainty [27] are shown as con- 
tinuous curves rather than error bars. 

The differences in C for every frequency are well 
within the bounds of even the best-case power meas- 
urement uncertainty, and show that the uncertainty due 
to the use of this interpolation method contributes only 
slightly to the overall uncertainty. 

We next examine the consistency of attenuation cor- 
rection coefficients determined from the two data sets 
1 1 years different in age. Figure 4 shows values of C 
based on values of a from all test results reported 
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Fig. 3. Attenuation correction coefficients for a from single-frequency 2005 data and least-squares fit of data 
for 24 frequencies. 
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Fig. 4. Attenuation correction coefficients for ideal and empirical values of a. 
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herein for 1994 and 2005, and also for the attenuation 
coefficient for ideal conditions [31]. For all frequen- 
cies, differences in C from the 1994 and 2005 data sets 
are small compared to power measurement uncertain- 
ties. This agreement is taken to confirm that the proce- 
dures for determining C are frilly adequate for the 
intended purpose. 

It is also clear that values of C based on empirical 
and ideal values of a differ insignificantly for frequen- 
cies below 10 MHz, and substantially for frequencies 
above 15 MHz. To the extent that similar behavior 
applies to standard source output power measurements 
made by other RFBs equipped with conical targets, 
similar plots could assist the determination of the fre- 
quency range for which attenuation corrections are 
necessitated by constraints on overall uncertainty. 



5. Conclusions 

An expression has been derived for the in-situ atten- 
uation correction factor C for a 45° conical target used 
to measure the output of a collimated, circularly sym- 
metric ultrasound source. 

In-situ attenuation a for distilled water was deter- 
mined for 24 frequencies between 4.8 MHz and 29.2 
MHz by least-squares linear regression analysis of 
measurements made in 2005 of ultrasound power ver- 
sus source-target distance. This analysis also confirmed 
the validity of the assumption that a is independent of 
position in the volume of water defined by the ultra- 
sound beam between the transducer and the NIST RFB 
target. 

The combined accuracy of the derived expression 
and reproducibility of the performance of targets built 
to the NIST design was assessed by re-analysis of data 
from similar measurements made in 1994 with a differ- 
ent target. For the two data sets spanning 1 1 years, the 
average values of R differed by 0.0007, and had an 
average value of 0.9970, indicating good accuracy and 
reproducibility. 

The results of least-squares linear regression analysis 
of the frequency dependence of a were found to con- 
firm the expected strong correlation between a and/^, 
and therefore to validate an interpolation method for 
determining a for arbitrary frequencies between 4.8 
MHz and 29.2 MHz. 

For this frequency range, the attenuation correction 
factor C was computed for both in-situ and ideal values 
of the attenuation of water. The difference between the 
two was found to be insignificant for frequencies below 
10 MHz, and was found to exceed the best-case com- 



bined relative expanded power measurement uncertain- 
ty of NIST RFB power measurements [27] for frequen- 
cies greater than 1 5 MHz. 
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